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METHOD FOR LOCATING AN ELEMENT OF INTEREST CONTAINED 



BACKGROUND OF THE INVENTION 

The invention relates to the location of an element of interest 
contained in a three-dimensional object on the basis of the positions of 
homologous regions of interest corresponding to the said element of 
5 interest and appearing in a set of stereotaxic images of the said object. 

Although the invention is applicable to any set of stereotaxic 
images of an arbitrary three-dimensional object, it is particularly useful in 
the medical field and, in particular, during a stereotaxic examination in 
mammography in order to locate microcalcifications in a breast. 

10 A stereotaxic examination in general makes it possible to access, 

with high precision of the order of one millimeter, a given point of a 
three-dimensional object, on the basis of two two-dimensional projections 
of this object acquired, for example, at two opposite angular incidences on 
either side of the normal to the plane of the image receiver delivering the 

15 said projections on the basis of the knowledge, on the one hand, of the 
position of the projection of this point into the plane of each stereotaxic 
image obtained and, on the other hand, of the geometry of the stereotaxic 
imaging apparatus which has led to the construction of the two images, it 
is possible to calculate the exact spatial position of this point in the 

20 three-dimensional object, by trigonometric calculation. 



IN A THREE-DIMENSIONAL OBJECT 



In summary, a stereotaxic examination generally requires: 



at least two images of the object at different angulations, and 
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perfect knowledge of the acquisition geometry of the stereotaxic 
system, the ability to locate the site of the projection of the chosen 
element of interest on the various stereotaxic images obtained. 

The first two requirements do not pose a major problem. Currently, 
5 stereotaxic examination in mammography is carried out with the aid of a 
mammograph equipped with a stereotaxic imaging device. The 
mammograph comprises an X-ray tube located at the end of a first arm 
which is mobile about an axle and emitting X-radiation to a receiver 
located at the end of another arm. A breast support plate, or patient 
10 support plate, on the one hand and, on the other hand, a compression 
plate which holds the breast in place during the mammography, are 
arranged between the tube and the receiver. The image receiver may be 
a digital receiver, such as a CCD camera, for example, delivering digitized 
stereotaxic images. The stereotaxic imaging requires the X-ray tube to be 
15 rotated around the breast support and compression plates, in two 
successive opposite orientations on either side of its initial position 
perpendicular to the plane of the image receiver. 

On the basis of the digitized stereotaxic images, obtained, the 
location of the element of interest on the various images is a difficult 
20 operation, in particular in mammography, since it requires corresponding 
homologous regions of interest in two different stereotaxic images to be 
set in correspondence or matched with a single element of interest 
located in the breast. 

Thus, in mammography, an attempt is made to match the regions 
25 of interest which correspond to the same microcalcification on the various 
stereotaxic images obtained. 
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However, the matching of the projected microcalcifications is a 
difficult problem because the microcalcifications do not look the same 
from one image to another. Their shape and their contrast may change, 
as may their arrangement in space. They may be superposed with fibrous 
5 zones or with other microcalcifications. 

To date, the matching has been carried out entirely manually by the 
radiologist. The latter would chose those regions of interest which were 
correspondence on the two angulated stereotaxic images, then enter the 
two-dimensional coordinates of these regions into a computer which gave 
10 the spatial position of the microcalcification, that is to say the three- 
dimensional coordinates of the latter. On the basis of these coordinates, 
the radiologist could adjust, for example, the position of a needle holder 
so that the needle becomes aligned with the microcalcification, in 10 order 
for it to be punctured, for example. 

15 However, this entirely manual matching proves to be a time- 

consuming and sometimes inaccurate step. 

BRIEF SUMMARY OF THE INVENTION 

The purpose of the invention is to substantially reduce the duration 
of a stereotaxic examination and to obtain better matching precision. 

20 In an embodiment of the invention, a similarity measurement is 

defined, making it possible to evaluate the resemblance between a region 
of interest (microcalcification, for example) selected (manually by the 
radiologist, for example, or else automatically by an automatic selection 
method) on one of the angulated stereotaxic images and a region of 

25 interest of a second stereotaxic image. 
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On the basis of this similarity measurement, a number of 
"candidates" belonging to a search zone and having a high degree of 
resemblance to the initial microcalcification are selected in order to obtain 
the homologous region of interest. 

5 Next, by using the principles of stereotaxy, the three-dimensional 

coordinates of these candidates are calculated, then they are reprojected 
onto the centering image. The similarity between these reprojected points 
and the microcalcification initially chosen is again calculated. In fact, if the 
candidate actually corresponds to the initial microcalcification, then its 

10 calculated projection will actually correspond to the image of this 
microcalcification on the centering image. 

Finally, with the aid of these two similarity values calculated for 
each candidate, the best candidate is proposed. 

in theory, it would be possible to omit the step of reprojection using 
15 the centering image and to make do with the matching as has just been 
explained in order to select the region of interest homologous to the one 
selected on the target image. However, it is by far preferable to use the 
reprojection of the candidates onto the centering image in order to verify 
the first selection made, in particular in the field of mammography in which 
20 the microcalcifications may change shape and contrast from one image to 
another. 

In other words, the invention is proposes a method for locating an 
element of interest contained in a three-dimensional object on the basis of 
the positions of homologous regions of interest corresponding to the said 
25 element of interest and appearing in a set of stereotaxic images of the 
said object. This method includes a step of selecting, in a first stereotaxic 
image, or target image, a first region of interest, in particular a 
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microcalcification, referred to as the target region of interest. The method 
also includes a matching of the said first region with a second region of 
interest, homologous to the first and appearing in a second stereotaxic 
image. 

5 According to a general characteristic of the invention, the 

stereotaxic images being digitized, a target pixel in the said target region 
of interest is selected in the selection step. In the matching step, a target 
window of chosen dimensional characteristics and containing the said 
target region of interest is generated around the selected target pixel. A 

10 set of pixels is then determined in the second image, according to a 
predetermined selection criterion,, and a second window, of the same 
dimensional characteristics as the said target window, is generated 
around each selected pixel. A correlation processing is carried out 
between the grey-scale levels of the pixels in each second window and 

15 the grey-scale levels of the pixels in the target window, so as to obtain a 
correlation value for each second window. The region of interest 
homologous to the target region of interest is then identified on the basis 
of the analysis of the set of correlation values thus obtained. The risks of 
matching error between homologous regions of interest are thus 

20 minimized. 

According to one embodiment, the analysis of the correlation 
values obtained includes the selection of a certain number of correlation 
maxima or minima, the homologous region being selected from those 
whose associated correlation value is one of these correlation maxima or 
25 minima. 

In fact, depending on whether the image is "normal" or inverted (a 
"normal" image being defined as having grey-scale levels and contrasts 
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which are similar to those usually visible on film negatives), use will 
respectively be made of maxima or minima. 

It is also possible to provide a determination of the dynamic range 
of the maxima or minima selected, and the comparison of the obtained 
5 dynamic values with a threshold. 

The correlation processing used may include a centered or 
uncentered normalized correlation processing, or else a normalized 
difference processing. 

As a variant to the matching explained above, the invention also 
10 proposes a matching based on the treatment of predetermined 
characteristics of any region of interest, in particular physical 
characteristics such as extension, contrast and gradient, then a 
comparison between these lists of characteristics associated with various 
regions of interest. 

15 More precisely, according to a general characteristic of the 

invention, a target pixel in the said target region of interest is selected in 
the selection step and, in the matching step, a target window of chosen 
dimensional characteristics and containing the said target region of 
interest is generated around the selected target pixel. On the basis of the 

20 grey-scale level values of the pixels in the said target window, a first set of 
numerical values, respectively representing predetermined characteristics 
characterizing any region of interest (for example shape, contrast or 
gradient characteristics, etc.) is determined for the said target region of 
interest, and this first set of numerical values is stored. In the second 

25 image, according to a predetermined selection criterion, a set of pixels 
belonging respectively to so-called candidate regions of interest appearing 
in the second images is determined, and a second window, containing the 
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said corresponding candidate region of interest, is generated around each 
of these selected pixels. For each candidate region of interest, a second 
set of numerical values, respectively, is representing the said 
predetermined characteristics, is determined on the basis of the 
5 grey-scale level values of the pixels in the associated second window, and 
this second set of numerical values is stored. The region of interest 
homologous to the target region of interest is then identified on the basis 
of a comparison processing between the first set of numerical values and 
each second set of numerical values. The risks of matching error between 
10 homologous regions are thus minimized. 

The comparison processing may include a distance minimization 
processing between the two sets of numerical values. 

In the two matching variants which have just been explained, and in 
order to further minimize the calculation time an epipolar zone containing 

15 at least the epipolar segment relating to the target pixel is determined in 
the second stereotaxic image. The pixels in the epipolar zone then contain 
the selected pixels. In order to reduce the calculation time further, instead 
of performing matching on all the pixels in the determined epipolar zone, 
the pixels in the second image are advantageously selected from those 

20 having grey-scale level maxima or minima whose dynamic range is 
greater than a predetermined threshold. 

Furthermore, in order to overcome the image background 
problems, and in order to further minimize the risks of matching errors 
between two homologous regions of interest, it is particularly 
25 advantageous for the matching according to the invention to include a 
prior filtering of the stereotaxic images, for example of the so-called 
"top-hat transformation" type. 
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The use of the step of reprojecting the candidate elements of 
interest into the centering image is independent of the type of matching 
previously used. 

Thus, according to a general characteristic of the invention, the 
5 location method of the type explained above comprises a selection, in a 
first digitized stereotaxic image, of a first so-called target region of 
interest, a selection, in a second digitized stereotaxic is image, on the 
basis of a first automatic matching, of at least one second region of 
interest which may be homologous to the target region, the determination 
10 of the spatial position of a candidate pixel of a candidate element of 
interest corresponding to these two regions of interest, the determination 
of the position, in a third stereotaxic image, of a projected pixel 
corresponding to the projection into this third stereotaxic image of the 
candidate pixel, and a second automatic matching between the target 
15 region of interest and a vicinity of the projected pixel, defining a projected 
region of interest. 

In practice, a set of second regions which may be homologous to 
the target region is advantageously selected in the second stereotaxic 
image, on the basis of the said first automatic matching, and the spatial 

20 position of each candidate pixel corresponding to each pair of regions of 
interest which is formed by the target region and one of the second 
regions is determined. The position in the third image of each 
corresponding projected pixel is also determined and the second 
automatic matching between the target region and a vicinity of each 

25 projected pixel is carried out. 

It is possible to order the set of second regions obtained in 
decreasing order of their probability of being the region homologous to the 
target region. It is also possible to order the set of projected regions 
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obtained in decreasing order of their probability of being the region 
homologous to the target region. The region homologous to the target 
region is then advantageously selected as being the one whose product of 
the respective ranks in the two orderings is a minimum. 

5 The first automatic matching may be one of the two matchings 

which have been explained above, that is to say based on a correlation 
processing or on a comparison processing. 

Similarly, the second automatic matching may include a correlation 
processing as defined above, between the vicinity of each projected pixel 
10 and a target window containing the target region of interest. This second 
automatic matching may also be a comparison processing as defined 
above, between the target region of interest and each projected region of 
interest. 

From a hardware point of view, the invention can be regarded as a 
15 device for locating an element of interest contained in a three-dimensional 
object on the basis of the positions of homologous regions of interest 
corresponding to the said element of interest and appearing in a set of 
stereotaxic images of the said object. This device includes a stereotaxic 
imaging apparatus equipped with a support for the said object and an 
20 image receiver such as a CCD camera capable of providing digitized 
stereotaxic images. Means for displaying these images are also provided, 
for example a video screen, such as that of a microcomputer. Selection 
means, for example a "mouse" of a microcomputer, make it possible to 
move a cursor on the screen and to "click" this marker onto a target pixel 
25 of a region of interest, making it possible to select the target pixel. Finally, 
a microprocessor and associated memories are provided in order to 
incorporate, in software form, means making it possible to functionally 
implement all the steps of the methods which have been described above. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Other advantages and characteristics of the invention will emerge 
on examining entirely non-limiting embodiments, and the appended 
drawings, in which: 

5 Figure 1 is a schematic view of three stereotaxic images of a 

three-dimensional object from a stereotaxic imaging device; 

Figure 2 very schematically illustrates the notion of an epipolar 
segment; 

Figure 3 illustrates an epipolar zone chosen according to the 
10 invention in order to perform a matching; 

Figure 4 is a general flow chart of a method of an embodiment of 
the invention; 

Figure 5 more particularly illustrates a step of a variant of the 
method, using a correlation processing; and 

15 Figure 6 schematically illustrates, more particularly, a method of an 

embodiment of the invention using reprojection. onto the centering image. 

DETAILED DESCRIPTION OF THE INVENTION 

A stereotaxic examination, in particular in mammography, is 
composed of a series of three exposures of a three-dimensional object 1 
20 (Figure 1), for example a breast, resting on a support 3 and compressed 
by a compression plate 2, with the aid of an X-ray tube which respectively 
occupies three different positions 6, 7, and 8. In practice, one exposure is 
taken at an angle of 0° and two exposures are taken angulated at two 
equal and opposite angles, in practice ±15°. 
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A CCD receiver 4 thus makes it possible to obtain three digitized 
stereotaxic images 10, 11 and 12. According to a convention which is 
normally used, the image 10 is the right image, while the image 12 is the 
left image and the image 1 1 is the centering image. 

5 An element of interest 5, for example a microcalcification, 

contained in the three-dimensional object 1 provides regions of interests, 
respectively referenced 5d, 5c and 5g, on each of the images 10, 11 
and 12. The centering image is used in particular to check that the 
microcalcification will be accessible during the possible subsequent 

10 medical intervention. 

The three-dimensional coordinates of the microcalcification 5 can 
be obtained, in a way which is conventional and well-known to the person 
skilled in the art, by a simple trigonometric calculation on the basis of the 
knowledge of the two-dimensional coordinates of the two projections 5d 
15 and 5g, for example. These coordinates are expressed with respect to an 
origin manifested by the center of a cross 9 etched on the plate 
supporting the breast, and the center of this reference cross is identified 
on each image in order to geometrically initialize the system. 

A person skilled in the art knows that, by definition, the epipolar 
20 segment 50 (Figure 2) of a point 5d of an image 9 of angulation 91, on an 
image 12 of angulation 02, corresponds to the set of points in the second 
image 12 which may have the point 5d as antecedent. In other words, it is 
the conical projection, by the source of angulation 92, of the straight-line 
segment 52 joining the said point to the source of angulation 91 . 

25 In the case where the stereotaxic images have been obtained with 

angulations such that the positions of the X-ray tube lie in a plane parallel 
to the lines of the images, this epipolar segment 50 is a straight-line 
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segment which is found to be parallel to the lines of the image. However, 
if the thickness of the three-dimensional object is not known, it is not 
possible to determine the position of the ends of this epipolar segment. 
This epipolar segment then becomes an epipolar straight line which, in the 
5 particular acquisition geometry explained here, lies on the same image 
line 12 as the point of interest 5d on the image 9. 

Furthermore, when use is made of a detector which has 
geometrical distortions, the epipolar straight-line is converted into an 
epipolar curve 53 (Figure 3). It is possible to determine the coordinates of 
10 the set of points on the epipolar curve on the image 12 by a geometrical 
calibration of the CCD camera. It would then be possible to perform the 
matching according to the invention only on this set of points. 

If, however, for reasons of precision of the geometrical calibration, 
the coordinates of the points on the epipolar curve can only be 

15 determined approximately, the choice will then be made to perform the 
matching according to the invention on an epipolar strip 51 which contains 
the epipolar curve. In practice, an epipolar zone, or search zone, bounded 
by two parallel lines passing respectively through two points augmented, 
on either side of the uppermost and lowermost points of the epipolar 53 

20 by a predetermined margin, will be defined. A margin of 3 to 5 pixels will 
generally be opted for. 

A more detailed description will now be given, with reference more 
particularly to Figures 4 and 5, of two types of matching in embodiments 
of the invention. 

25 The first step consists in the selection (step 400) of a target pixel 

60d within a target region of interest 5d present, for example, in the right 
image 10, referred to as the target image. In practice, this selection can 
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be made manually by clicking on a point of the region of interest in the 
image in question appearing on a screen of a microcomputer. This being 
the case, it is also possible to provide an automatic pixel selection 
method. This pixel 60d may or may not be centered in the target region. 

5 Next (step 401), a target window 61 d containing the said target 

region of interest is generated. This target window can be generated 
automatically after having applied a conventional region-increase software 
tool to the target pixel 60d. As a variant, depending on the field of 
application of the invention, provision may also be made to choose 

10 windows of predetermined size. In all cases, a rectangular or square 
target window will generally be obtained, the sides of which have 
dimensions of the order of ten to a few tens of pixels. 

The grey-scale levels of each of the pixels in this target window are 
then stored in a memory associated with a microcomputer which 
15 incorporates as software all the functional tools for implementing the 
described method. 

On this basis, a set of pixels in a second stereotaxic image, for 
example the left image, will be selected in a step 402. This set of selected 
pixels 'will constitute a search zone 51 within which one or more candidate 
20 regions of interest which may be homologous to the target region of 
interest will be selected. 

As explained above, and with the aim of reducing the matching 
time, this search zone 51 will, for example, be the epipolar zone 51 
described above. More precisely, a second window 61 g of the same 
25 dimensional characteristic as the target window will be generated (step 
403) around each of the pixels selected in the search zone. The 
grey-scale levels of the pixels in each second window are stored and a 
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correlation processing is carried out (step 404) between the grey-scale 
levels of the pixels in each second window 61 g and the grey-scale levels 
of the pixels in the target window 61 d. 

The correlation processing then makes it possible to perform a 
5 similarity measurement between the target window and each of the 
second windows of the second image, so as to obtain a similarity value 
(correlation value between these two windows). 

It is on the basis of the analysis of these correlation values that the 
target region of interest will be matched with its homologous region of 
10 interest. 

As regards to the correlation processing (step 404), use may be 
made of a normalized correlation. More precisely, each correlation value 
of a second window will then be given by the formula: 

N 

E^yw > 

K=l _ ~~~T^ # J 

15 

where N is the - number of pixels in a window, k the running index of a 
pixel, l(k) the grey-scale level of the pixel k in the target window and J(k) 
the grey-scale level of the pixel k in the second window. 

It is also possible to use a normalized difference, making it possible 
20 to obtain the correlation value by the formula: 
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ln theory, the homologous region of interest will be the one 
contained in the second window which has the maximum correlation 
value. This being the case, several correlation values may in practice be 
found with a high intensity. A restricted set of correlation maxima are then 
5 selected (step 405) from the set of correlation values, according to a 
predetermined criterion such as the intensity of the maximum or such as 
the number of maxima in the list of correlation values ranked in 
decreasing order of intensity. 

It is also possible to use, as an additional selection criterion, a 
10 determination of the dynamic range of the maxima selected and the 
comparison of the dynamic values obtained with a predetermined 
threshold. The determination of the dynamic range of maxima is perfectly 
known to the person skilled in the art. However, he may optionally refer to 
French Patent Application No. 91 15308 for more details regarding such a 
15 dynamic-range determination. It should be noted here that the method in 
the embodiments of the invention makes it possible to determine, and 
display on the screen of the microcomputer, for example by flashing, the 
region of interest considered as being the region homologous to the target 
region, and to do this with the minimum risk of error. 

20 This being the case, the method in the embodiments of the 

invention will in practice be able to select a restricted number of candidate 
regions of interest which may be homologous to the target region. The 
radiologist will then have the possibility of manually selecting the one 
which he considers to be the actual homologous region. 

25 As explained above, a specific comparison processing may be 

used instead of a correlation processing in the matching processing. 
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Thus, as more precisely illustrated in Figure 4, a list of 
characteristics allowing general characterization of a region of interest is 
determined (step 409) for the target region contained inside the target 
window 61 d, and on the basis of the grey-scale levels of the pixels 
5 contained in this target window. These characteristics may include, for 
example, a value representing the extension of the region of interest, a 
value representing its mean width, a value representing its mean 
grey-scale level, or its mean intensity gradient. 

The same operation (step 410) is performed for the candidate 
10 regions of interest contained in second windows 61 g defined on the basis 
of pixels selected in the search zone. In this regard, in order to limit the 
processing time, the pixels selected for the determination of the second 
windows encompassing the candidate regions are chosen from the pixels 
of the search zone which have grey-scale maxima whose dynamic range 
1 5 is greater than a predetermined threshold. 

On the basis of these various lists of characteristics thus obtained, 
a comparison processing 41 1 is carried out between the list of 
characteristics assigned to the target region and each of the lists of 
characteristics assigned to each of the candidate regions. 

20 In practice, although another type of comparison processing may 

be envisaged, use will preferably be made of a comparison processing 
which provides for the minimization of a distance, for example the 
Euclidean distance, between the list of characteristics of the target region, 
then considered as a vector, and a list of characteristics of a candidate 

25 region, also considered as a vector. 

By way of example, provision can also be made to determine the 
mean value of the absolute values of the differences between the 
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homologous values of two lists of characteristics, and to search for the 
lowest of the mean values in order to deduce the homologous region of is 
interest therefrom. 

The selection of the homologous region provided in step 405 is 
5 thus obtained, with this variant, by minimizing a distance between lists of 
characteristics, rather than by maximizing a normalized correlation or a 
normalized difference. 

The invention provides some further improvements. 

Thus, in the matching variant providing for a correlation processing, 
10 it is possible to limit the search zone to those pixels of the epipolar zone 
which have grey-scale level maxima whose dynamic range is greater than 
a predetermined threshold. 

It is furthermore preferable, in particular in order to overcome the 
image background problems, to filter the stereotaxic images beforehand, 
15 with filtering of the "top-hat transformation" type. This type of filtering is 
known to the person skilled in the art and for more details he may refer to 
the article by J. SERRA Image Analysis and Mathematical Morphology, 
Vol. 2, Academic Press 1988. 

This type of filtering makes it possible to extract the light and 
20 narrow regions of a digital image and makes it possible to overcome the 
background. 

The steps of the method in the embodiments of the invention which 
will now be described, and providing in particular for reprojection onto the 
centering image, thus, in particular, equate to a verification of the results 
25 of the preceding matching between the target image and the second 
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image, and thus make it possible to further minimize the risks of matching 
errors. 

This being the case, although the rest of the method, which will 
now be described with more particular reference to Figures 4 and 6, 
5 provides the possibility of using correlation or comparison processings of 
the type already described, this method according to the invention, using 
reprojection onto the centering image, is independent of the prior 
matching used between the target image and the second image. 

It is now assumed (Figure 6) that the method according to the 
10 invention has made it possible to select, on the left image 12, two 
candidate regions of interest 5g and 20g which may be the region 
homologous to the target region 5d appearing on the target image 10. 

On the basis of the target pixel in the target region of interest and 
of the candidate pixel which has given rise to the selection of the 

15 candidate region of interest 5g, a stereotaxic trigonometric calculation and 
knowledge of the geometry of the stereotaxic device are used to 
determine (step 406) the three-dimensional coordinates of the 
corresponding point of the candidate object of interest corresponding to 
these two regions 5d and 5g. Then, by inverse stereotaxic trigonometric 

20 calculation, the coordinates in the centering image of the projected pixel 
corresponding to the projection of the pixel of the candidate element into 
this centering image are determined. 

It can, therefore, be seen that, in the present case, since the 
candidate region 5g is effectively the region homologous to the region 5d, 
25 the position of the pixel projected into the centering image 11 actually 
corresponds to a projected region of interest 5p which is actually the 
projection of the element of interest 5 into this centering image. 
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In contrast, using the same calculations, the candidate region of 
interest 20g will lead to a projected position 20p in the centering image 
being obtained which, in the present case, will not be superposed with any 
region of interest physically present in this centering image. 

5 Thus, a correlation processing between the target window 

containing the target region 5d and a vicinity of the projected pixel, 
containing the projected region of interest 5p, will give a correlation value 
greater than the correlation value obtained between the target window 5d 
and a vicinity of the projected pixel 20p. The same would have been true 

10 even if the vicinity of the projected point 20p had corresponded to a 
projected microcalcification other than the one actually corresponding to 
the microcalcification 5. 

This method therefore makes it possible to select, from the two 
candidates 5g and 20g, the region 5g which is actually the region 
15 homologous to the region 5d. 

Of course, the conclusion would have been the same when using a 
comparison processing 412 instead of a correlation processing 407. 

The use of a vicinity of the projected pixel, typically a few pixels 
around the projected point, makes it possible to account for the possible 
20 errors due to a shift in the object between the exposure at 0° and the 
exposures at ±1 5°. 

More generally, the set of candidate regions obtained in the second 
image can be ordered in decreasing order of their probability of being the 
region homologous to the target region. Thus, the region which has, for 
25 example, given the correlation maximum with the greatest dynamic range 
will be assigned rank 1 , and so on in decreasing order. The same ordering 
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can be performed for the set of projected regions obtained on the 
centering image. The region homologous to the target region will then be 
selected as being the one for which the product of the respective ranks in 
the two orderings is a minimum. 



5 



Various modifications in structure and/or function and/or steps may 
be made by one skilled in the art without departing from the scope and 
extent of the invention. 
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WHAT IS CLAIMED IS: 

1 . A method for locating an element of interest contained in a 
three-dimensional object comprising the steps of: 

a. providing a set of digitized stereotaxic images having 
5 positions of homologous regions of interest corresponding to the element 

of interest and appearing in the set of stereotaxic images; 

b. selecting in a first stereotaxic image a first target 
region of interest having a target pixel; 

c. matching the first region of interest with a second 
10 region of interest homologous to the first region of interest and appearing 

in a second region of interest; 

d. matching a generated target window of chosen 
dimensional characteristics and containing the target region of interest 
around the selected target pixels; 

15 e. determining a set pixels in the second image 

according to a predetermined selection criterion so as to generate a 
second window having the same dimensional characteristics as the target 
window around each selected pixel; 

f. processing a correlation between the gray-scale 
20 levels of the pixels in each second window to obtain a correlation for each 

second window; and 

g. using the correlation values to identify the region of 
interest homologous to the target region of interest and :thereby minimize 
the risk of matching error between the homologous regions of interest. 
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2. The method of claim 1 wherein the correlation values 
includes the selection of a certain number of correlation maxima or 
minima, the homologous region of interest being selected from those for 
which the associated correlation value is one of the correlation maxima or 
minima. 

3. The method of claim 2, wherein the analysis of the 
correlation values obtained includes a determination of the dynamic range 
of the maxima or minima selected, and the comparison of the dynamic 
values obtained with a threshold. 

4. The method of claim 1, wherein the correlation processing 
includes a normalized correlation processing. 

5. The method of claim 1 wherein the correlation processing 
includes a normalized difference processing. 

6. The method of claim 1 wherein the matching includes a prior 
filtering of the stereotaxic images. 

7. A method for locating an element of interest contained in a 
three-dimensional object comprising the steps of: 

a. providing a set of digitized stereotaxic images having 
positions of homologous regions of interest corresponding to the element 
of interest and appearing in the set of stereotaxic images; 

b. selecting in a first stereotaxic image a first target 
region of interest having a target pixel; 

c. matching the first region of interest with a second 
region of interest homologous to the first region of interest and appearing 
in a second region of interest; 
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d. matching a generated target window of chosen 
dimensional characteristics and containing the target region of interest 
around the selected target pixels; 

e. determining and storing a first set of numerical 
5 values, respectively representing predetermined characteristics of any 

region of interest for the target region of interest on the basis of the 
gray-scale values of the pixels in the target window; 

f. determining in the second image a set of pixels 
belonging respectively to candidate regions of interest appearing in the 

10 second image according to predetermined selection criterion; 

g. generating a second window, containing the 
corresponding candidate region of interest around each of the selected 
pixels; 

h. determining and storing a second set of numerical 
15 values, respectively representing the predetermined characteristics for 

each candidate region of interest, on the basis of the gray-scale values of 
the pixels in the associated second window; and 

i. identifying the region of interest homologous to the 
target region of interest on the basis of a comparison processing between 

20 the first set of numerical values and each second set of numerical values 
so as to minimize the risk of matching error between the homologous 
regions of interest. 

8. The method of claim 7, wherein the predetermined 
characteristics comprise shape, contrast or gradient characteristics. 
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9. The method of claim 7 wherein the comparison processing 
includes a distance minimization processing between the two sets of 
numerical values. 

10. The method of claim 8 wherein the comparison processing 
5 includes a distance minimization processing between the two sets of 

numerical values. 

1 1 . The method of claim 1 wherein an epipolar zone containing 
at least the epipolar segment relating to the target pixel is determined in 
the second stereotaxic image and the pixels in the epipolar zone contain 

10 the selected pixels. 

12. The method of claim 7 wherein an epipolar zone containing 
at least the epipolar segment relating to the target pixel is determined in 
the second stereotaxic image and the pixels in the epipolar zone contain 
the selected pixels. 

15 13. The method of claim 1, wherein the selected pixels are 

chosen from the pixels in the second image which have a gray-scale level 
maxima or minima whose dynamic range is greater than a predetermined 
threshold. 

14. The method of claim 7, wherein the selected pixels are 
20 chosen from the pixels in the second image which have a gray-scale level 

maxima or minima whose dynamic range is greater than a predetermined 
threshold. 

15. The method according to claim 7, wherein the matching 
includes a prior filtering of the stereotaxic images. 



14XZ00014 



16. A method for locating an element of interest contained in a 
three-dimensional object comprising the steps of: 

a. providing a set of digitized stereotaxic images having 
positions of homologous regions of interest corresponding to the element 

5 of interest and appearing in the set of stereotaxic images; 

b. selecting in a first stereotaxic image of a first target 
region of interest; 

c. selecting in a second stereotaxic image on the basis 
of a first automatic matching of at least one second region of interest 

10 which may be homologous to the target region; 

d. determining the spatial position of a candidate pixel of 
a candidate element of interest corresponding to the two regions of 
interest; 

e. determining in a third stereotaxic image of a projected 
15 pixel corresponding to the projection into the third stereotaxic image of the 

candidate pixel; 

f. providing a second automatic matching between the 
target region of interest and a vicinity of the projected pixel; and 

g. defining a projected region of interest so as to 
20 minimize the risk of matching errors between the homologous regions of 

interest. 

17. The method of claim of claim 16 wherein the selecting in the 
second stereotaxic image is on the basis of the first automatic matching a 
set of second regions which may be homologous to the target region; 

25 determining the spatial position of each candidate pixel 

corresponding to each pair of regions of interest which is formed by the 
target region and one of the second regions; 

determining the spatial position in the third image of each 
corresponding projected pixel; and 
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performing the second automatic matching between the 
target region and a vicinity of each projected pixel. 

18. The method of claim 17 wherein the set of second regions 
obtained is ordered in the decreasing order of their probability of being in 

5 the region homologous to the target region; 

wherein the set of projected regions obtained is ordered in 
the decreasing order of their probability of being the region of homologous 
to the target region; and 

wherein the region homologous to the target region is 
10 selected as being the one whose product of the respective ranks of the 
two orderings is a minimum. 

19. The method of claim 16, wherein the first automatic 
matching is a correlation processing matching. 

20. The method of claim 16, wherein the second automatic 
15 matching includes a correlation processing matching between the vicinity 

of each projected pixel and a target window containing the target region of 
interest. 

21. The method of claim 16 wherein the second automatic 
matching includes a comparison matching between the target region of 

20 interest and each projected region of interest. 
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METHOD FOR LOCATING AN ELEMENT OF INTEREST CONTAINED 
IN A THREE-DIMENSIONAL OBJECT 

ABSTRACT OF THE DISCLOSURE 

The stereotaxic images being digitized, a target pixel in a target 
5 region of interest is selected, a target window of chosen dimensional 
characteristics and containing the said target region of interest is 
generated around the selected target pixel, a set of pixels is determined in 
a second image, according to a predetermined selection criterion, a 
second window, of the same dimensional characteristics as the said target 
10 window, is generated around each selected pixel, a correlation processing 
between the grey-scale levels of the pixels of each second window and 
the grey-scale levels of the pixels of the target window is carried out so as 
to obtain a correlation value for each second window, and the region of 
interest homologous to the target region of interest is identified on the 
15 basis of the analysis of the set of correlation values thus obtained, so as 
to minimize the risks of matching error between the homologous regions 
of interest. The element of interest is then located on the basis of the 
positions of the two homologous regions. 
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